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1. Introduction 

The scaling dimension A for twist two operators at large values of the Lorentz spin S 
is characterized by the universal scaling function (cusp anomalous dimension) f{g): 

A-S = fig) \nS + ..., 



with g = gyM^ f^"^- The logarithmic dependence of the dimension on large 
Lorentz spin is a generic feature that has been independently observed for both 
M = 4 SYM and the corresponding string theory dual, see, e.g., |l|-H]. Importantly, 
such a scaling behavior stems 0, ^ from the large spin limit of the Bethe equa- 
tions [@]-0 which underlie the integrable structures of gauge and string theories. 
A linear integral equation arising in the large spin limit from the postulated set of 
gauge/string Bethe equations has been recently derived in ||10[. We will refer to it 



as the BES equation. This equation allows one to compute the universal scaling 
function f{g) to any desired order in perturbation theory. For instance, for the first 
few orders one obtains 

m = 8g' - ^ttV + . . . 
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Remarkably, independent gauge theory computations 12| confirm this resuh up 
to four loop orders, i.e. up to g^. 



Schematically, the BES equation has the following form 

a{t) = l + g^ j K{t,t')a{t')dt' , 



where K{t,t') = K^"'\t,t') + KK^''\t,t') and n = 2. Here is the so-called 

main scattering kernel and K^^'^ is the dressing kernel. The dressing kernel arises 
upon taking into account an additional scattering phase in the S-matrix which is 
apparently needed to achieve an agreement with string theory predictions in the 
strong coupling limit This phase is also the only freedom in the S-matrix which 
cannot be fixed by symmetry arguments |jTB|.^ The form of the phase is largely 
but not completely constrained by an additional requirement of crossing symmetry 



. The constraints arising from crossing symmetry were shown to hold |jT8[ at two 



leading orders in string perturbation theory [0, [T^. Recently an interesting solution 



for the crossing symmetric phase has been obtained pO| which agrees nicely with all 
available string and field-theoretic data. 

The agreement observed at the four loop level between the scaling function of 
BES and the results obtained in the field-theoretic framework is very important, 
because it implies the appearance of the additional scattering phase in perturbation 
theory and, therefore, it distinguishes the BES equation from the earlier proposal 
which corresponds to taking k = 0. 



According to the AdS/CFT duality |21] the large g behavior of f{g) can be 



extracted by studying the energy of the folded string spinning in the AdSa part of 
the target space (the GKP solution) and it turns out to be 0, |^ 

31n2 , , 

f{g)=^g + ... (1.1) 

TT 



In the recent work and an important question was addressed, namely, how the 
perturbative expansion controlled by the BES equation can be extrapolated to large 
values of the coupling constant and wether the resulting expression is consistent with 
the string theory predictions. In particular, in |^ a numerical method was devised 
to solve the BES equation around the strong coupling point and numerical solutions 
were shown to be in perfect agreement with eq. (|1.1|) , providing a non-trivial test of 
the AdS/CFT correspondence. 

In this note we find an exact analytic expression for the density o"(t) in the strong 
coupling limit. If one insists on truncating the BES equation at leading order for 
large g the kernel becomes degenerate and therefore a unique solution for o"(t) cannot 
be found without additional assumptions. However, for finite values of g the solution 



^See p4|-[p6[ for recent studies of the string S-matrix and its symmetries. 
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is unique and can be found numerically with high precision. These two statements 
are consistent because expanding the BES equation in a power series in l/^f we find 
a second equation at subleading order that removes any ambiguity in the strong 
coupling solution ait). 

To gain some intuition our approach will be to first investigate the solution nu- 
merically for finite values of the coupling, which will suggest a simple additional re- 
quirement that should be imposed on the underdetermined, truncated BES equation 
to single out a unique solution in the strong coupling limit which is the convergence 
point of our numerics. After this additional requirement is identified, we find an ana- 
lytic solution for a{t) in the strong coupling limit. We then show that this solution is 
in fact completely determined by taking into account the subleading contributions to 
the BES equation, thus confirming analytically the validity of the auxiliary condition 
suggested by the numerics and of our strong coupling solution cr(t). 

Upon making the Fourier transform to the rapidity u-plane the corresponding 
density a{u) reads 



\u\-r 



V2 



\ 



1 



1-^ / 



where 6{u) is the step function. Thus, on the rapidity plane the leading density 
is an algebraic function which is constant in the interval |m| < 1. We see that, in 



contrast to the weak coupling solution of the BES equation [|T2], the strong coupling 
density exhibits a gap between [— oo, — 1] and [1,cxd], see Fig.l. Remarkably, this 
is reminiscent of the behavior of the corresponding solutions describing classical 
spinning strings. 



-0.5 



Fig.l. The plot of the analytic solution for the leading density Ang^aiu). 



We also would like to point out that the strong coupling solution of the BES equation 
we have found is apparently different from the singular solution Xsing discussed in 

ra. 
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The paper is organized as follows. In section 2 we analyze the matrix form of the 
BES equation at strong coupling, first numerically and then analytically. We compute 
the coefficients in the expansion of a{t) in terms of Bessel functions and based on 
this numerical analysis we make a guess of how these coefficients could be (partially) 
related to each other at strong coupling. In the strong coupling limit we obtain an 
analytic equation for the coefficients which turns out to have a degenerate kernel. 
Supplementing this equation with the proposed relation among the coefficients allows 
us to find an exact analytical solution for cr{t). We then show that in fact no guess 
is necessary if one expands the BES equation to higher order in 1/g, which leads to 
a second equation that singles out our solution as the correct one. 

We proceed to find an analogous pair of equations for the subleading coefficients 
at strong coupling. Again they appear at different orders in the inverse coupling 
constant, each of them individually being a degenerate, half-rank equation. We 
identify some constraints on the subleading value of crit), but postpone a detailed 
analysis of the strong coupling perturbation theory to a future publication. In section 
3 we analyze the BES equation at strong coupling in the Fourier space and verify 
that the results of section 2 agree with this analysis. Finally, we end with conclusions 
and open problems. 



2. BES matrix equation at strong coupling 

In this section we study the density of fluctuations ait) for large values of the coupling 
constant. First we argue, by performing numerical analysis, that a{t) obeys a certain 
additional requirement. Then we consider the matrix BES equation at leading order 
in the large g expansion, use the previously found requirement to solve for the density 
and conflrm by expanding the BES equation further that this solution is in fact 
complete determined analytically. Finally, we derive and briefly discuss constraints 
on the subleading value of cr{t). 

2.1 Numerical studies of a(t) 

The universal scaling function is related to a{t) as follows 

f{g) = lQgMO), (2.1) 
where the density of fluctuations a(t) satisfles the following integral equation |10| 



a{t) = {K{2gt, 0) - Ag^ dt' K{2gt, 2gt') a(t')) • (2.2) 

The kernel separates into two pieces, the main scattering kernel K^™) and the dressing 
kernel K^^^ 

Kit, t') = K^"^^ [t, t') + 2ir(") (t, t') , (2.3) 
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which are expressed in terms of Bessel functions as follows 

_ Jiit)Mt')-Jo{t)J,{t') ^ 

poo J.// 

K^'^\t, t') = Ag' / dt" K,{t, 2gt")^^K,{2gt\ t' 
.In e — i 



Here Kq and Ki denote even and odd parts of the kernel under change of the sign 
of the arguments: 



oo 



Ko{t, t ) = ^2 _ f2 = ^ Z^(2n - l)J2n^l{t)J2n-l[t ) 



n=l 
oo 



^ ej,flj„m-tMt)J.<f ) ^ I ^ 2,a.«)..„(0 . (2.4) 



n=l 



As was shown in \22], if one introduces a function 



s{t) = ^-j^a{t) (2.5) 
and expands s{t) in the following series 

sit) = ±s.^ (2.6) 

n=l ^ 

then eq.( |2.2| ) transforms into a matrix equation for the vector s = (si, S2, S3, ...)* 

(l + ir(™) + 2ir(^))s = (l + 2C)e. (2.7) 

Here e = (1, 0, 0, 0, ...)* and K^"'\ K^'^l C are infinite-dimensional matrices with the 
entries 

K^Z^ = 2(iVZ)^„ , 

Cmn = 2{PNZQN)mn, 

where Q = diag(l, 0, 1, 0, ...), P = diag(0, 1, 0, 1, ...) and N = diag(l, 2, 3, ...). The 
coupling constant enters into the matrix Z whose entries are given by 

Z„„.rd«^«^. (2.8) 

Jo He* - 1) 

The key observation is that for intermediate values of g, for instance 2 < g < 20, 
one can approximate the infinite-dimensional matrices entering the BES equation by 
matrices of finite rank d, with d not much larger than g. With matrices of finite 
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rank it is possible to solve numerically for the coefficients Sk for different values of 
the coupling constant and to find the best fit result for an expansion of the type 



Sk 



1 



1 

? 



(2.9) 



As the numerical analysis indicates, the finite rank approximation is valid for comput- 
ing the coefficients Sk with k <^ d. We have solved numerically eg. ( ^.71 ) for numerous 
points in the range 2 < g < 20 and use d = 50. In the table below the values for a 
few leading coefficients are exhibited. 



k 


S2k-1 


'^2k 


1001(4-1 


1 


0.500006 


0.499993 


0.003 


2 


-0.75005 


-0.74977 


0.038 


3 


0.93727 


0.93676 


0.055 


4 


-1.09281 


-1.09415 


0.12 


5 


1.2239 


1.2333 


0.77 



Some comments are in order. First, notice that the value for s{ is in perfect agreement 
with the value predicted from string theory s{ = 1/2 and confirmed numerically 
(with a precision higher than the one presented here) by |2^. Second, notice that 
the difference between s^^-i and is in all the cases smaller than 1%, and gets 
bigger as k increases; this is related to the fact that the rank of the matrices is finite. 

Thus, the numerical analysis suggests that in the limit of infinite rank matri- 
ces the following relation holds for the leading coefficients in the strong coupling 
expansion 

4-1 = 4- (2.10) 

As we will see later on, this condition will allow one to solve analytically for the 
coefficients and the values obtained will be in perfect agreement with the ones 
computed numerically. 

Further evidence comes from the fact that one can approximate the matrix ele- 
ments Zmn by their analytic values at strong coupling (see next subsection). There- 
fore, one can consider matrices of much higher rank, fix a sufficiently large value of 
g and compute (numerically) the coefficients s^. Below we present the results for 
g = 10000 and d = 250. 



k 


J 

S2k-1 


S2k 


1001(4-1 -4)/4l 


1 


0.49993 


0.49991 


0.0049 


2 


-0.74943 


-0.74938 


0.0073 


3 


0.93585 


0.93577 


0.0085 


4 


-1.09033 


-1.09023 


0.0089 


5 


1.22455 


1.22444 


0.0087 
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As d increases, we see that the difference between S2i^_i and (for k = 2,3,4) 
decreases considerably. Also, the difference is approximately constant for small values 
of k. We should stress, however, that a priory there is no reason to expect the results 
obtained by keeping in the BES equation only the leading term for Z to be valid, 
since as we will see the subleading terms in the matrix elements Zmn are necessary to 
fix uniquely the leading order solution for Sk- Surprisingly, one still obtains a good 
approximation to the large g solution in this fashion. Nevertheless, we regard the 
present computation as less robust. 

To conclude, requiring continuity in g, the leading coefficients si in the strong 
coupling expansion of the function a{t) exhibit the relation (|2.10|) , which constraints 
a possible form of 

2.2 Analytic solution at strong coupling 

For finite real values of g the matrix element Zmn is given by a convergent integral. 
However, it is not obvious if it is possible to express the result of integration as a 
power series in 1/g. Indeed, expanding the integrand in eq. as 







t \t 2 2Ag 



leads to a power series (with coefficients given essentially by the Bernoulli numbers) 
which converges only for \t/ {2g) \ < 2tt. We see that only the first two leading terms in 
this expansion can be integrated, while already in the third term divergent integrals 
appear. It is therefore natural to assume that the expansion of Z develops as 

Z = gZ^ + Z'^ + (2.11) 

where the first two terms, Z^ and Z^^, are given by the convergent integrals mentioned 
above. Explicitly, they are 

« 8 cos((m — ?T,)7r/2) 



TT (m + n + l)(m + — l)(m — n + l)(m — n — 1) ' 
_ 1 sin((m-r2)7r/2) 
TT — n'^ 

Assumption (|2.11|) is well supported by the numerics. We have checked that numeri- 
cal values for Z^^ for large g are in a very good agreement with the analytic expres- 
sions ( p.l2| ). On the other hand, the status of the higher order terms in eq.( p.llj ) is 



not clear to us^. Thus, in what follows, we largely restrict our investigation of the 
BES equation to the first two leading terms. 



^For instance, whether expansion ( 2.11 ) contains terms of the type \og{g)/g, etc. From the string 
theory point of view, the appearance of such terms could be an indication of a two-loop divergence 
of the world-sheet S-matrix. Further analysis is needed to clarify this issue. 



-7- 



With this word of caution we proceed to investigate eq. ( |2.7|) in the strong 
couphng hmit. It is not hard to see that it imphes the following equation for the 
leading vector 

K^-)^s^ = C^e, (2.13) 

where the leading matrices K^"^^ ^ and are obtained by keeping the leading contri- 
bution only. It turns out that for even values of d the kernel K^'^^^ has rank d/2, 
hence from the equation above it is possible to solve only for half of the components 
of s. However, now the condition ( |2.10| ) enters into play. Imposing eq. (|2.10| ) together 
with eq.( p.l3| ) allows one to uniquely determine the vector s^. 

In order to find the solution in the limit of infinite d it is convenient to express 
the leading equation in the following way 

K°s° + K^s^ = ^e, (2.14) 

where s°, are vectors of length d/2 comprising the odd and even components of 
respectively: s° = {s{, s|, Sg, ...)* and = {s^, si, Sg, ...)*. Further e = (1, 0, 0, ...)*. 
It is easy to check that eqs.( |2.14] ) and ( |2.13| ) are equivalent provided 

(K )mn = -Z^2m-l,2n-l + ^2m+l,2n~l 7 (-^ )mn = ^2m-l,2n + ^2m+l,2n ■ (2.15) 



Then, we find that the unique solution of eq. (|2.14|) satisfying the relation s" = 
turns out to be 

4..^4^(-if«^|i|. (2.16) 

This remarkably simple expression for the coefficients is the main result of this 
paper. By using the following identities (true in the limit of infinite rank) 

{K^ ~^)^n = -4(-l)™+"+Vr^^ n<m 

= -4(-l)™+"+^m3, n>m (2.17) 
32 

^ ^ (4^2_(i_2^)2)(i_2n)2 

one can check explicitly that the coefficients (p.l6|) indeed solve eq. (|2.14|) . Thus, 
we found that restricting the coefficients Sk to the leading order expressions the 
function s{t) is given by 

where we have omitted the subleading contributions. The last formula can be con- 
sidered as the leading term in the large g expansion of the density s{t, g) with gt kept 
finite. As is clear from eq.(2^), we are only interested in values of s{t) for t > 0. 
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The series can be summed and for this range of t the result expressed in terms of the 
confluent hypergeometric function of the second kind U{a,b,x): 

s{t) = (r(|)f/(-i,o,-4^^?t) + r(f)f/(i,o,-4z^t)) 

+ ^^e-'''\n^U{-l^A^9t) + T{^U{\,^Atgt)). (2.19) 

The leading density ( |2.19|) has a rather complicated profile. For (jft ^ it perfectly 
reproduces the desired result s{t) 1/4(7. On the other hand, the asymptotic of 
s{t) for gt —>■ oo exhibits a highly oscillating behavior. 

2.3 An alternative derivation of the strong coupling solution 



Let us now show that the result ( p.l6| ) can in fact be derived without resorting to any 
auxiliary conditions obtained from numerical arguments. The degenerate equation 
appearing at leading order, which can be used to express one half of the coefficients 
Sn in terms of the other, can be supplemented by another half-rank equation from 
subleading terms in the BES equation. Together they determine a unique solution. 

Writing out the integral equation (|2.2| ) in the basis (|2.6| ) with explicit matrix 
indices we find 

Jn{2gt) _ Ji{2gt) , ^ ^ J2n{2gt) Jn{2gt) 

— 1" "'^^2?i,l ^nZjnmSm ;:r~ 

2gt 2gt 2gt 2gt 

-16n(2m - l)^2n,2m-1^2m-l,rSr^^^^^ , (2.20) 

where all indices are summed over from 1 to oo. Now we split up the integral equation 
according to powers of g and into odd and even rows (indices of Bessel functions). 

At 0{g) the odd equation is trivial and the even one reads: 

2(2m — 1)Z2„ 2m-1^2m-l,r'Sr = ^2n,l = ■ (2.21) 

This is precisely eq. (|2.13|) employed in the previous subsection. At C(l) the odd 
rows lead to the condition ^ 

^2m-l,r^i = 2^m,l ■ (2.22) 

Actually this equation implies the previous one. It determines one half of the coef- 
ficients in terms of the other. Expanding further, at 0{1) the even equation is 
given by 

8'^^2n,l ~ 4nZ2„ „jS^ — 16n(2m — l)^2n,2m-1^2m-^l,r'5r^ 

- 16n(2m - l)ZL,2™-i^2l-i,.4 - 16n(2m - l)Zf^,2m-iZL-i,rsi = • (2.23) 
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To determine the other half of the coefficients we need to ehminate from this 
equation. To do this we examine the odd equation at 0{l/g) 

-2(2m - - 2(2m - l)Zf^_,,y, = sL_, . (2.24) 

Now we use (|2.22| ) and ( p.24| ) to simphfy (|2.23|) , which gives 



^2n,m^m '^^2n,2m-1^2m-l ~ ^2n,m( 1) — ^ ■ (2.25) 



This is the second equation we were looking for, which together with ( |2.22 ) com- 



pletely determines s^. If we define a matrix Znm which is identical to Znm except 
for a sign fiip when both n is even and m is odd, we can combine the two conditions 
into the full-rank equation for the strong coupling solution 

ZLsL = ■ (2.26) 

Note that the additional minus signs in the definition of Znm arise precisely because 
of the introduction of the dressing kernel, and would by absent for k = 0. Writing 
out ( ^.261 ) explicitly shows that the have to satisfy 



E 



E 



4(-l)"+^+i4,„, 


(2n 


-2k- 


l){2n 


-2k + l){2n + 2k- 


3)(2n 


+ 2k- l)7r 








^2n 












8n{2n - 1) 


8{n - 


l)(2n- 1) 














{2n 


-2k- 


l){2n 


-2k + l){2n + 2k- 


l){2n 


+ 2A; + l)7r 



-5n,i , (2.27) 



W + "2n-l .2.28) 



8n(2n + 1) 8n{2n 



for all n > 1 (where it is understood that the second term on the right hand side of 
( p.27|) is absent for n = 1). Indeed these equations are obeyed by 



-s^ _ (zl)!:H2n-l)!! 

S2n-1-S2n- 2"(n-l)! ' ^ ' 

which is precisely the solution (|2.16|) found in the previous subsection. To show this, 
note that the coefficients S2„_i are generated by the Taylor expansion of (1 + x)~'^/^, 
which makes it easy to perform the above sums as integrals over functions of the 
form x™(l — x^)^^/'^ for some appropriate power m chosen to generate the necessary 
terms in the denominators of the left hand sides of eqs.( 2.27 ) and ( |2.28 ). 
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2.4 Fluctuation density in the rapidity plane 

To get more insight into the structure of the leading solution, we find it convenient 
to perform the (inverse) Fourier transform of the density (j{t) a{u): 



p / dt e'2f*«e-l*l/V(|t|). (2.30) 



We recall that u is a rapidity variable originally used to parameterize the Bethe root 
distributions of gauge and string theory Bethe ansatze p, which gives another 
reason for studying the density of fiuctuations on the w-plane. Thus, substituting in 
eq. (|2.30|) the power series expansion for a{t) we get 



1 _ _ 

— Ysn dt e^^gt. -|t|/2. 

n=l 

For large values of g this expression can be well approximated as 

^ ^ ^^ti J^oo 2g\t\ 
The last integral is computed by using the following formula 

Jo 2gt n \ \ J 

with = 1/2 =1= 2igu. 

In this way we obtain the following series representation for the density (T(n) 

^ oo 
^ n=l 

where 

/n = ^ [K(1 + Vl + V/M'))"" + K(l + Vl + V/(u-)2))-'^" . 

In what follows it is convenient to introduce the expansion parameter e = 1/ (2g). 
Our considerations above suggest that the density cr{u) expands starting from the 
second order in e: 

a{u) = ^a\u) + e-a'\u) + . . . (2.31) 

To find the leading contribution o^[u) we have to develop the large g expansion of 
the functions /„. The result is not uniform, it depends on whether n is even or odd 
and also on the value of u . For n even we find 

^T2,,(M) + 0(e) for |m| < 1 , 



2k 



2k 

2k 



(-1)' 
2k 



uh+Jl-^)) +0{e) for \u\>l 
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Here T2k{u) are the Chebyshev polynomials of the first kind. For n odd we obtain 

i^v^r^ U2k-2{u) + C(e) for |m| < 1 , 



2k- 



2k-l 



+ C(e) for |n| > 1 



Here U2k-2{u) are the Chebyshev polynomials of the second kind. We recall that 
the Chebyshev polynomials of the first and the second kind form a sequence of 
orthogonal polynomials on the interval [—1, 1] with the weights (1 — v?)"^^"^ and 
{1 — u^Y^"^ respectively. 

To find the leading density a^{u) inside the interval |n| < 1 we can use the 
trigonometric definition of the Chebyshev polynomials which corresponds to choosing 
parametrization u = cos 6 with < 6 < n. Thus, taking the limit g ^ oo we obtain 
for the leading density the following expression 



cr [u 



k=l 



(2.32) 



where 



f2k{9) = (-1) 



.cos2k6 
2k 



f2k-i{9) = -{ 



_fcSin(2A^ 



2A; - 1 



(2.33) 



Given the proposal (|2.16|) for the coefficients s^, we can now sum the series for a^{u) 
and obtain 



1 



a [u) = -. 

TT 



(2.34) 



Thus, inside the interval [—1,1] the density a^{u) is constant. 



-3 -2 



-□.25 
-0.5 
-□.75 



Fig. 2. The left figure is the plot of a numerical solution for the leading density 
TTCT^ (u) . The right figure represents an exact analytic solution for the same quantity. 



We note that the functions fk{G) provide a complete set of functions on the interval 
[0, 2tt], while on the interval [0, vr] they are overcomplete. This fact is consistent with 
the ambiguity for the coefficients Sk we observed in the previous section. Moreover, 
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it is not hard to convince ourselves that the proposal for discussed in the previous 
section is the only possible solution such that Sj^.i = and for which the density 
a^{u) is constant for \u\ < 1. 

Further, it is easy to sum up the series defining the leading density for \u\ > 1. 
The result is 

a\u) = 7^ I 2 - I . (2.35) 



271 



1 - u{u =F Vu"^ - 1) 



Here the minus and plus signs in the denominator corresponds to the regions u > 1 
and u < —1 respectively. The plot of the complete analytic solution for a^{u) is 
presented in Fig. 2. It should be compared to the plot of the numerical solution 
obtained by using in the series representation for the density the coefficients 
obtained from our numerical analysis. The numerical plot corresponds to taking 
g = 4:0 and truncating the series at k = 20. 

Finally, we mention the expression for the scaling function f{g) in terms of the 
density a{u): 



oo 



f{g) = 32g' / dua{u) = 8g / dua'{u) + ... = % + ... 

J —oo J —oo 

This completes our discussion of the leading order analytic solution of the BES 
equation in the strong coupling limit. 

2.5 Subleading corrections 

Here we will investigate the first subleading correction to the leading coefficients s\.. 
By expanding the BES equation we have already obtained equation ( p.24| ) which 
expresses the subleading density s*^ in terms of the leading one: 

^2rn-l,r^r — ~^2rn-l,r-^r " 2(2m — 1) "^2™"^ ' (2.36) 

Again this equation is degenerate and needs to be supplemented by a second one 
that appears at higher order in \ j g in the BES equation. Substituting the explicit 



form of the coefficients s„ obtained above, the right hand side of eq. ( |2.36| ) is seen to 
vanish (i.e. = 0), which already implies tight restrictions on the subleading 

corrections s'^^ . As discussed above, this equation allows for solving for half of the 
components. For instance, using eqs.( p.l4"|) and (|2.17| ) we can solve for the even 
components in terms of the odd ones 

32 

= - g(-l)'"-"- ^^^2 _ (1 _ 2^)2)(i _ 2n)2'2"-i • (2.37) 
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To obtain another equation constraining the subleading solution, we examine 
corrections to the BES equation as follows. From the 0{l/g) contribution to the 
even rows of the BES equation ( p.20| ) we find 



'2n 



(2.38) 



16n(2m - 1) 



yl ss£ I yl ysl si _i_ yl yssl 

1^2m-l,r^r ^2n,2m-1^2m-l,r^r ^2n,2m-1^2m-l,r^ 



2n,2m- 



7SSI J 



I ysi yi si _i_ yst ysl i , ysst yl £ 

"t~^2n,2m-1^2m-l,r''r ^2n,2m-1^2m-l,r"'r ^2n,2m-l'^2m-l,r^r 



To eliminate the term in s***^ we turn to the odd rows of the 0{l/g^) equation 



-2(2m - 1) 



ssl I yst St I y 



■'2m- 



sst t 
2m—l,rr 



'2m-l 



(2.39) 



Using this together with ( |2.22| ) and ( |2.24| ) to simplify ( |2.38| ) we obtain 



^2n,m\ 



-irst 



_ys£ (T^m I ^ £ 

^2n,m\ ^) "m 4^ 



(2.40) 



This can be combined with (p.36|) into the full rank equation 



7« s£ 

nm^m 



_7se e ^ £ 

nni^m 2?7. ' 



(2.41) 



where again the right hand side vanishes when evaluated on the leading solution 
found above. 

In view of the divergence of the 0{l/g) correction to Znm, which could indi- 
cate the presence in it of terms e.g. 0{log{g)/ g), the strong coupling expansion of 
Sn{g) may contain more than just pure powers of 1/g. This complication begins at 
0{l/g^) in Sn and suggests that the equations for the subleading terms have to be 
supplemented by additional constraints. Luckily, the equations for the leading term 
in section 2.3 are not affected by this complication. We hope to analyze the structure 
of subleading terms more fully in future work. 



3. Leading density from the integral equation 

In the appendix we compute the inverse Fourier transform of the BES equation 
(rather than that of its leading solution) to the rapidity variable u.^ In this section 
we will expand the BES integral equation in the small parameter e = 1/ {2g) assuming 
that the density expands as in eq.( p.31[ ) and rederive the leading order solution a^{u) 



by different means. 

It turns out that the branch cut of the square root functions defining the all- 
loops Bethe-ansatz forces us to distinguish the regimes \u\ ^ 1. The leading order 



■^This approach to the BES equation was also pursued by A. BeUtsky (seminar at Princeton 
University, December 2006). 
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BES equation constrains the leading density cr^(u) only within the unit interval; no 
condition is obtained for \u\ > 1. We will show that this equation has a unique 
solution, namely the constant function a^{u) = l/vr, \u\ < 1. 

In the discussion of the last section, the even and the odd functions, f2k and 
/2fc-i, constitute equivalent complete sets of functions on the interval |m| < 1. The 
constraint that the density be constant in this region allows one to eliminate, e.g., 
the coefficients S2fc_i in favour of (or vice versa). Now, to leading order only the 
even functions f2k have support outside the interval \u\ < 1. The knowledge of the 
density outside this interval would then fix the remaining coefficients S2j.- 

Above, we have rather taken the opposite approach: pairing the numerical ob- 
servation that S2fc_i = with the requirement that the density be constant within 
the unit interval we could equivalently fix the density cr^(-u) on the whole real axis, 
c.f. eq.(g:35D. 

The leading terms in the strong coupling expansion of the BES equation are 
obtained in appendix 5.2. We find 



1 



(3.1) 



where 



u, u 



'1- u 



i\2 



( 



V 



+ ^log 



(u'y 



V 



w 



(3.2) 



In Fourier space, the potential for the BES equation is actually the kernel with the 
second argument put to zero; or to put it differently, the kernel integrated on a delta 
function. In the rapidity plane this means that the potential is the integral of the 
kernel on a constant function. By construction, our leading equation is solved by 

1 



a [u] 



\u\ < 1 



TT 



(3.3) 



Let us check whether there is a second solution inside the interval \u\ < 1. To this 
end, we put 



u 



V 



1 



'\2 



and differentiate the whole equation in U. Next, we substitute f/^ 
to obtain 



1 + 



1 

2U 



dva{v) 



1 

u 



u-v 



(3.4) 
U and -^V 

(3.5) 



where a is the density written in the new variables but including the transformation 
of the integration measure: 

a = a , = (3.6) 
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There is no regular density that could yield 1/U from the Hilbcrt transform (i.e. the 
second term on the r.h.s.). We will therefore look for a density that produces the 
constant term on the l.h.s. from the Hilbert transform, and whose norm is defined 
by the 1/U terms. Now, 

and in the first case the normalization is 1/2, while in the second it is —1/2. Note 
that the symmetric combination of the two trial densities 

^ = (3.8) 



goes to zero under the Hilbert transform whereas it does affect the norm. Thus, a 
solution exists and it is unique. We have to choose 

1 \/V 

a{V) = - / (3.9) 

which corresponds to cr{u') = l/vr. 

In appendix 5.2. we further derive the next-to-leading order of the BES equation. 
We find 



\u\ <1. 



By what was said above the leading dressing kernel is invertible on [—1,1]. We 
therefore conclude 

a'\u) = 0, |m| < 1. (3.10) 

No constraint is found on the subleading density outside the unit interval. 

We refrain from pushing the analysis any further for the following reasons: First, 
at the next order we would perhaps encounter the same degeneracy while we lack an 
independent condition like s^i^ = Sg^.i that would allow us to make progress. Second, 
the naive expansion in e quickly leads to expressions with rather hard singularities 
in u, u' which are difficult to handle consistently. 

The original BES equation has the property that the energy is given by the value 
of the density at zero. In the rapidity variables this means that we can recover the 
leading contribution to the energy from the normalization of a^. Alternatively, we 
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may use the formula 



E{9) 



log(5)8^2 
21og(5) 



1-2^2 / dua{u) 



x^{u) x~{u) 
1-^^ dua{u) (2VI -u^-e + ..^ - 

- — y dua{u) 



+ 



It is instructive to see how the value of the energy predicted from string theory is 
reproduced: The contribution of the constant part of within the unit interval 
cancels the leading 1 in the square bracket. On the other hand, the next-to-leading 
term would receive a contribution from an a priory non- vanishing subleading density 
inside the unit interval. We get 



E{g) = log{S) Sg 



du 



- 2 dua'\u)VT^^ 



+ . . . = log(5) 4g + ... 



Here we used an identity 



/ dua[u)— j du j=J= = / dua[u). 

Jo Ji u^(^l + ^1- ^) -1 J--^ 

satisfied by the leading density and also the fact that the subleading density van- 
ishes inside the unit interval. In fact, the absence of the subleading density can be 
considered as the consistency test on our equations. 



4. Conclusions 

In this paper we analyze the strong coupling limit of the BES equation which de- 
scribes the universal scaling function of high spin operators in A/" = 4 gauge theory. 
We have shown that expanding the BES equation in inverse powers of the coupling 
constant leads to two equations for the large g solution, one appearing at leading the 
other at subleading order in 1/g. Together they determine a unique solution in the 
strong coupling limit whose exact analytic form we present. 

Obviously, the next step would be to understand the structure of the higher order 
perturbation theory around the strong coupling point and, in particular, to derive 
the next-to- leading corrections to the universal scaling function f{g). We make some 



-17- 



progress in this direction by deriving a pair of equations for the subleading solution, 
but also point out some subtleties that arise in developing the expansion around the 
leading order solution. 

As was shown in ||T0| (see also |^), the coefficients of the perturbative series 



describing the solution of the BES equation at weak coupling admit an analytic 
continuation to strong coupling, where they coincide with those predicted by string 
theory. The approach we adopt here can be considered as another, complementary 
way to analytically continue from weak to strong coupling. 

On the rapidity w-plane the leading fluctuation density a^{u) appears to be 
constant inside the unit interval \u\ < 1. We could argue that this constant part of 
a^{u) is an artifact of the way the BES equation was derived: The non- vanishing 
constant part of the leading density offsets the splitting of the weak-coupling density 
into a log divergent one-loop part and a regular higher loop piece carrying log 5* 
as a coefficient. Further, as we have shown, the subleading correction inside the 
unit interval is absent; in a manner of speaking a gap opens between [—00, —1] and 
[1,cxd]. This could be qualitatively compared to the results obtained from string 
theory. Indeed, the solution of the integral equation describing the classical spinning 
strings in AdSs x in the limit S 00 with spin J along fixed has support 

only outside the interval \u\ < 1. The same behavior is expected^ for the GKP 
solution which is obtained in the limit J ^ 0. For finite S the solution is elliptic 
and it exhibits logarithmic singularities in the limit S 00. On the other hand, 
our strong coupling density ( |2.35| ) is an algebraic function which carries log 5* as a 
normalization. Of course, this density leads to the same energy as for the GKP 
string. Thus, it is desirable to understand the detailed matching between the string 
density (higher conserved charges) and the density we found from the strong coupling 
limit of the BES equation. We plan to return to this interesting question elsewhere. 



Acknowledgements 

We are grateful to Andrei Belitsky, Sergey Frolov and Matthias Staudacher for valu- 
able discussions. I. R. K. and M. K. B. would like to thank Sergio Benvenuti and 
Antonello Scardicchio for collaboration on the early stages of this project. This work 
was supported in part by the EU-RTN network Constituents, Fundamental Forces 
and Symmetries of the Universe (MRTN-CT-2004-005104), by the INTAS contract 
03-51-6346, by the NWO grant 047017015, and by the National Science Foundation 
under grant No. PHY-0243680. The work of L. F. A. was supported by the VENI 
grant 680-47-113. The work of G. A. was supported in part by the RFBI grant 
N05-01-00758 and by the grant NSh-672.2006.1. 

^We would like to thank Sergey Frolov for the discussion of this point. 



- 18 - 



5. Appendix 

5.1 The inverse Fourier transform of the BES equation 

Let us define the kernels 

1 + 



1 + ^ 



X (u)x+{u') I \ X (u) X (u') 



1 + .-..L.... 1 - 



, X (u) x'^ (u') I \ X (u) X (u') , 

K.{u,u') = ^ \ i ) , ( , (5.2) 



with 



^Hu) = ^{u±-] |l + 1 • (5.3) 



We further define 

Ko{u,u') = '-du\og(^K+{u,u')K_{u,u')') , (5.4) 

Kiiu,u') = '-d^\og(^K4u,u') / K4u,u')^ . (5.5) 
In the same way as in appendix D of |^ we may show that 

ko,i{u,u') = 2/ / / dtdtV"*+»'*'|t|e-(l*l+l*'l)/2iro,i(2^7|t|,2^7|t'|), (5.6) 



where Kq^i are defined in {2A) in the main text. The dressing kernel is |10] 

dt"iri(2(7|t|,2(7|t"|)^;^^o(2<?|t"|,2<?|t'|) (5.7) 
and its Fourier back-transform yields 

f poo 

K^''\u,u') = 2g^ / / dtdtV"*+'"'*>|e-(l*l+l*'l)/2ir(")(2^l^l,2^l^1) 

J J —OO 

= ^J J J dvdv'dt"ko{u,v)j-^k,{v',u'). (5.8) 
In the u, u' variables the BES equation becomes: 

= 2na{u) - 2 / du' aiu)- 

+ J du' ((T(n') - ^) (i^o(M, u') + (m, u') + 2i^(^) (m, u'] 



(5.9) 



We have written the potential as the integral of the kernels on the constant function 
— l/(47r5f^), i.e. the Fourier back-transform of —1/(2(7^) 5{t'). 
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5.2 The strong coupling limit 

The integration over t" in eg. ( ^.81) is not well defined. To make sense of it we borrow 
the derivative from Kq (c.f. ( |5.4|) ) on the right: 



^ dt"^--^ = dt"e-(-^')*" + t" j (5.10) 

= (-^(i(i; + v')) - ^(-2(t; + v')) + 2ti5{v + i;')) , 

where \&(w) is the digamma function. We may then partially integrate the outer 
derivative back to act on Kq. 

The physical situation we consider is in a kinematical regime in which u scales 
with 2g at strong coupling, as can be seen for instance from the numerical studies 
of the integral equation for density of fiuctuations in In the following we will 
rescale u —>■ 2gu and consider a perturbation series in e = l/{2g). 

In the middle integral in eq.( |5.8| ) we rescale and employ the asymptotic expansion 
of the digamma function: 

-*(i(t; + v')) - <iJ{-i{v + v')) + 2tt5{v + v') (5.11) 

^ - log ({v + v'Y) +-Siv + v') + 0{e^) 

9 

Note that we have discarded a term log(25f) + cq, where cq is an integration constant 
from the differentiation trick. A constant integrated onto the leading order of Ki in 
the second argument (or Kq in the first) actually yields zero. Further, the derivative 
of the exact expression in terms of the digamma function is a principal value for 
l/{y + v') (i.e. a cut-off at zero). We usually evaluate log((f + v'Y) under the 
subsequent f , v' integrations by partial integration, albeit using a Cauchy principal 
value. It ought not matter which definition is used for the principal value. 

After rescaling we find the following leading order dressing kernel in configuration 
space: 

k^^)\u,u') = -^l 1^ dvdv'kiiu,v) 01og((t; + t;')'))^oK,«')- (5-12) 
Here i^o,i = e-^o,i- To leading order 




\u\>l : x+{u) = x-{u) = —u + Y 1 - ^ j + • • • (5-13) 

|m| < 1 : a;+(u) = ^^{u + - u^^ + . . . (5.14) 
^ ^-u - - m2 j + . . . 



X [u] 
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from which it is easy to obtain 



K^{u,v) = : |t;|>l, Ki{u,v) = : |u| > 1 (5.15) 
and the non- vanishing sectors 

|u| < 1 : Ko{u, v) = n{S{u - v) + 5{u + v)) , 



(5.16) 



|m| > 1 : Kq{u,v) 



1 - 4 



— 9 9 ' 

1 — 



and 



l-u] < 1 : K^{u, v) = tt{6{u - v) + 6{u + v)) - 



(5.17) 



\v\>\ : K[{u,v) 



1_1 



y2 y2 _ y2 



In particular, the leading contribution ( |5.12|) to K^^\u,u') is only non- vanishing 
when both u,u' lie in the unit interval. Let us evaluate eq.( 5.12| ) explicitly, plugging 
in the formulae for Kqi and executing the two integrations in turn. First, 



— 1 I'OO 

+ / 



dv'hog{{v + v'f) k'oW,u') 



= -2tt log (^v^l - {u'y + Vl-v^ 
= -TV log {v^ - (u'y) : 

while the f '-integration over the unit interval contributes 



: \u'\, \v\ < 1 
n'l < 1 , If I > 1 , 



(5.18) 



Hence, if Im'I, If I < 1 



oo 1 

df'-log((f + f')') K',iv',u') 



(5.19) 



(5.20) 



-2tt log (Vl - {Wy + VT^) + \ log ((f 2 - (m')2)2) 



TT 



log 



1 - [u 



'\2 



and otherwise we find zero. It is then easy to see that 



\ 



log 



'i-(«o')'y 



'1 - (u 
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to leading order in the strong coupling expansion. 
The rescaled BES equation is 

= 2Tia{u) - 2 / du'aiu')- (5.21) 

Note that the kernel in the second term on the r.h.s. is a representation of vr — m') 
so that the first two terms cancel for small e. The numerical analysis displayed in 
the paper suggests to look for a solution 

The leading equation is 

du' [o\v!) - i:^^) ^(u, u') = . (5.22) 

It constrains the density only inside the unit interval. We prove in the main text 
that the obvious solution 

a\vi) = - : < 1 (5.23) 

TT 

is also the only one. At the next order we find the equation 

+ f\u'a'\u')k^'^\u,u'). (5.24) 



Now, to leading order Kq{u,u') vanishes when \u'\ > 1 so that the Kq term can be 
dropped. We show below that the other two kernels in the first line cancel. Hence 
for any the last equation reduces to 

1 
1 



It follows that 



: \u'\ < 1 (5.25) 



because K^*^^^^ is invertible on the unit interval (see Section |^). At this stage we 
cannot make a statement about cr*^ outside the unit interval. 

To show the aforementioned cancellation of the two kernels, let us work out 
u,u') for \u'\ > 1. In this case, Kq[v',u') in (|5.12| ) vanishes, so that we pick 
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up the subleading correction in the right factor while the other two terms may still 
be taken at leading order. From eq. (|5.4|) we can work out 



\u\ < 1, \u'\ > 1 : u') = , (5.26) 

\u\ > 1, \u'\ > 1 : K'o\u,u') = ^ , 
where we have introduced the notation 



A = 1)(-2 + m' + K)') 



It is not too hard to check that 

r>00 



oo 



\v\ < 1, \u'\ > 1, (5.27) 



= : \v\,\u'\>l 
from which we may easily deduce that 

i^W^^(M,M') = _iir[(M,M') : \u\ < 1, \u'\ > 1. (5.28) 
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